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ABSTRACT 

Context. Multiwavelength (MWL) observations of the blazar PKS 2155-304 during two weeks in July and August 2006, the period when two 
exceptional flares at very high energies (VHE, E> 100 GeV) occurred, provide a detailed picture of the evolution of its emission. The complete 
data set from this campaign is presented, including observations in VHE y-rays (H.E.S.S.), X-rays (RXTE, CHANDRA, SWIFT XRT), optical 
(SWIFT UVOT, Bronberg, Watcher, ROTSE), and in the radio band (NRT, HartRAO, ATCA). Optical and radio light curves from 2004 to 2008 
are compared to the available VHE data from this period, to put the 2006 campaign into the context of the long-term evolution of the source. 1 
Aims. The data set offers a close view of the evolution of the source on different time scales and yields new insights into the properties of the 
emission process. The predictions of synchrotron self-Compton (SSC) scenarios are compared to the MWL data, with the aim of describing the 
dominant features in the data down to the hour time scale. 

Methods. The spectral variability in the X-ray and VHE bands is explored and correlations between the integral fluxes at different wavelengths 
are evaluated. SSC modelling is used to interpret the general trends of the varying spectral energy distribution. 

Results. The X-ray and VHE y-ray emission are correlated during the observed high state of the source, but show no direct connection with longer 
wavelengths. The long-term flux evolution in the optical and radio bands is found to be correlated and shows that the source reaches a high state 
at long wavelengths after the occurrence of the VHE flares. Spectral hardening is seen in the SWIFT XRT data. 

Conclusions. The nightly averaged high-energy spectra of the non-flaring nights can be reproduced by a stationary one-zone SSC model, with only 
small variations in the parameters. The spectral and flux evolution in the high-energy band during the night of the second VHE flare is modelled 
with multi-zone SSC models, which can provide relatively simple interpretations for the hour time-scale evolution of the high-energy emission, 
even for such a complex data set. For the first time in this type of source, a clear indication is found for a relation between high activity at high 
energies and a long-term increase in the low frequency fluxes. 

Key words. Galaxies: active — BL Lacertae objects: individual: PKS 2155-304 — Radiation mechanisms: non-thermal — Gamma rays: galaxies 
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1. Introduction 

Several aspects of the exceptionally high state of the blazar 
PKS 2155-304 that was observed by the Imaging Air Cherenkov 
Telescopes (IACT) of the High Energy Stereoscopic System 
(H.E.S.S.) d uring the summer of 200 6 ha ve already been 
presente d in lAbramow ski et al. I d2010l) and I Aharonian et al. I 
( 2009b]). These publications focused on the spectral and tem- 
poral variability in the very high energy (VHE, E>100 GeV) 
band and on the short-term multiwavelength (MWL) behaviour 
during the night of the second VHE flare (hereafter "Flare 2"; 
MJD 53946). The present paper completes this work by provid- 
ing a long-term MWL view from radio to VHE during the whole 
of the 2006 campaign and beyond. The MWL behaviour of the 
source during the 2006 high state is analysed with stationary and 
time-dependent emission models and put in the context of the 
long-term data set. 

Obtaining a multispectral view of the temporal evolution 
of blazar fluxes is one of the main keys towards a global 
understanding of AGN physics. MWL observations help to 
gain detailed insight into the acceleration and emission pro- 
cesses of relativistic particles, which are thought to occur 
close to the central black hole ( e .g. iNeronov & Aharonian I 
l2007t IRieger & Aharonian I 12001 ISol & Istomin I 12009k 
in the jet (e.g. ISikora et al. 1 119941: llnoue & Takaharal 



Katarzvn skt Sol & Kus I 1200 ll 120031: iTavecchio & Ghisellini 



W) or 
1996; 
sellini 



2008). The double-bump structure of the spectral energy distri- 



bution observed in blazars at high energies is interpreted differ- 
ently in leptonic and hadronic models, which suppose that the 
particle population in the source is dominated by either electrons 
and positrons or by hadrons. 

The first bump, which is found at X-ray/UV energies in high- 
frequency peaked BL Lacs (HBLs) and at optical energies in 
low-frequency peaked BL Lacs (LBLs) and flat-spectrum radio 
quasars (FSRQs), is usually ascribed to synchrotron emission 
from a relativistic population of electrons. The second bump, at 
GeV/TeV energies in HBLs and usually at hard X-rays in LBLs 
and FSRQs, is interpreted as inverse Compton (IC) emission 
from electrons up-scattering synchrotron or external photons in 
leptonic models. In the hadronic scenarios, it arises instead from 
proton synchrotron emission or from hadronic interactions in- 
side the jet or between the jet outflow and the ambient medium. 

For HBLs, external photon fields are generally assumed to 
play a minor role in the emission, since no strong emission lines 
are observed. Thus the synchrotron self-Compton (SSC) sce- 
nario, where the bump at highest energies is attributed to syn- 
chrotron photons that have been up-scattered by their parent 
electron population, is often applied to interpret spectral energy 
distributions from these objects. SSC models require only a very 
restrained set of free parameters, compared to those leptonic sce- 
narios that add an external component or compared to the more 
complex hadronic models. Throughout the modelling presented 
in this paper, the SSC approach will be adopted and the term 
"one-zone model" will be applied to models where both the X- 
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ray and y-ray components are largely dominated by emission 
from the same zone. 

If one wants to put constraints on the different existing sce- 
narios, the active states of blazars provide crucial information on 
the characteristics of the emission region and the emission pro- 
cess, since they can exhibit high variability on short time scales 
and over a wide wavelength range. 

PKS 2155-3 04 is one of the brightest BL Lacs in 
the X-ray (e.g. iBrinkmann et al. I 11994 iKubo et al. I Il998t 
iGiommi et al~l Il998h IVestrand & Sreekumarl 1 19991) and EUV 
(Marshall, Carone & Fruscione 119931) bands. The source is 
classified as an HBL and has a redshift of z — 0.116. 
Since the discovery of X-ray emission from this object 
(Griffi ths etalTI 1 19791 ISchwartz et all Il979h . it has been re- 
peatedly observed over a wide range of frequ e ncies from ra- 
dio to VHE y-rays (e.g. iTreves et al. I |1989t lEdelsonetal 
1 1995b IVestrand. Stacv & Sreekumarjl 19951 IZhang & Xie lll996t 
Urrvetal.l 11997) iPian etal I 119971; iPiner & Edwards I 



2004 



Dominici et al. I 12004c I Aharonia n et al. I I2005at IZhang et al. I 



2006; Osterm an et al. I [2007). BeppoSAX observations from 
1996 to 1999 show X-r ay variability on time scales of « 1 
hour (Zhang et al. 2002). Several physical impl ications for the 



emission mechanisms have been reported (e.g. Chiarjoetti et al. 
1999j: |Kataoka et al. Il2000t lEdelson et al. 1120011 iTanihata et al. 
boOlclzhang et al. 1120021 120061) . 

In November 1997, an outburst from PKS 2155 
304 was detected in X-ray s and y-ray s with EGRET 
(ISreekumar & Vestrandl I1997D. BeppoSAX jC hiappetti et al. I 



1999). and RXTE ( IVestrand & Sreekumarl [1999J). During this 
active phase, the Durham group reported the first detection of 
VHE y-rays at a level of 6.8 standard deviations (cr) above 300 
GeV dChadwick et al. Ill999allbh . In 1998, when the X-ray flux 
level was low, they had no t found any evidence of TeV y-ray 
emission ( Chadwick 1 1 1 9991) . PKS 2155-304 was also observed 
with the CANGAROO-I 3.8m tele scope in 1997 No y-ray 
signal above 1.5 TeV was detected dRoberts et al. 1 1X999). The 
source was further observed in 1999, 2000, and 2001 with the 
CANGAROO-II telescope. It remained in a low state of X-ray 
activity in those periods, and was not detected above the energy 
threshold of 4 20 GeV dNishiiima et al. 11200 ll iNishiiima | [2002: 
lNakasell2003h . 

PKS 2155-304 was confirmed as a TeV y-ray source by the 
H.E.S.S. collaboration with observations in 2002 and 2003. A 
detection with a sig nificance of 45 cr at ener gies greater than 160 
GeV was reported ( Ah aronian et al. Il2005al) . A first MWL cam- 
paign of PKS 2155-304 including H.E.S.S. data was conducted 
in 20 03 over several weeks, du ring the construction phase of the 
array (I Aharonian et al. f2005bl) . The source was observed simul- 
taneously in the X-ray range by RXTEfPCA, in the optical by 
ROTSE and in the radio band by the Nan^ay Radio Telescope 
(NRT) and was found to be in a low state. Intra-night variabil- 
ity was seen in the VHE and X-ray band, with the shortest time 
scales detected with RXTE «25 minutes. The optical flux showed 
only moderate variation. No correlation was observed between 
the different bands. The hardness ratio in X-rays showed an in- 
crease in the spectral hardness with higher flux levels. Data from 
another MWL campaign conducted by the H.E.S.S. collabora- 
tion in 2004 were affected by the poor quality of the atmosphere 
at the time of data acquisition. 

The most recent MWL observations of the source in 2008, 
including data from H.E.S.S., ATOM, RXTE and for the first 
time data i n the high-energy y-ray r ange from Fermi, have been 
reported in lAharonian et al.l d2009ah . The source was found with 
a relatively low VHE flux level and no correlation between VHE 
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y-rays and X-rays was observed. There was, however, some ev- 
idence of correlated behaviour between the VHE and optical 
bands. 

The observation of a spectacular flare (hereafter "Flare 1") 
on July 28, 2006 (MJD 53944), with a more than 20-fold flux 
increase in one night a nd variability time scale s down to 200 
seconds was reported in Aharonian et al. (2007). Observations 
with Cangaroo-III confirmed the detecti on of a very active stat e 
with variability on sub-hour time scales ( Sakam oto et al. 1 12008). 
The very fast variability seen by H.E.S.S. implies in the one- 
zone SSC framework that the Doppler facto£] dh had a value of 
about 60 to 120 times the size of the emission region in units 
of the Schwarzschild r adius of the central black hole (see also 
Begelman et al. 1 (120081) 1. This means that either the Doppler fac- 
tor was unusually large compared to previous blazar observa- 
tions or that the emission region was much smaller than the size 
of the central black hole. 

A study of simultan eous observations in the X -ray and VHE 
bands during Flare 2 (Aha ronian et al. ] |2009b) found a very 
steep correlation between the X-ray flux detected by CHANDRA 
and the VHE y-ray flux, with the flux in the VHE band de- 
creasing approximately as the third power of the X-ray flux. 
No indication for a time-lag between those bands was detected. 
As a consequence of the steep correlation, one-zone SSC mod- 
els are strongly disfavoured to describe the emission in X-rays 
and VHE y-rays during this flare (for a general discussion see 
iKatarzvnski et al. 1120051) . 

A further investigation of the emission mechanism in the 
VHE band during the four most active nights, MJD 53944 to 
53947, has led to the conclusion that the obse rved flux variability 
stems from a lognormal stochastic process dAbramowski et alTI 
2010). The existence of a quiescent state of the source and 
of complex spectral variability in the VHE band was equally 
demonstrated by combining the data from 2006 with older and 
newer data sets. 

In the following, the complete two-week MWL campaign 
of 2006 is described and an interpretation of the observed flux 
evolution in the different energy bands is given. The evolution 
of the VHE flux and s pectrum of PKS 2155-30 4 on different 
time scales, covered by Abramow ski et al. I d2010l) . is here com- 
pared to the MWL emission from the radio to the X-ray band. 
The physical p rocesses behind Flare 2 were discussed in gen- 
eral terms by Aharonian et al. (2009b). Here, the data from this 
event are put into the context of the MWL observations before 
and after the flare and are confronted with actual emission mod- 
els. A comprehensive summary of all the available MWL data 
that correspond to the period of the H.E.S.S. campaign is pro- 
vided in Section|2] In Section|3] a study of correlations between 
light curves at different wavelengths is presented and the spectral 
variability in the X-ray and VHE band is described. Modelling 
of the nightly averaged spectral energy distributions with a sta- 
tionary one-zone SSC model is presented in Section [4] For the 
night of Flare 2, the only night with simultaneous H.E.S.S. and 
CHANDRA coverage, three different time-dependent SSC mod- 
els are discussed. The insight gained from these models and from 
the MWL analysis is discussed in Section|5] No attempt is made 
at describing the very rapid variability observed in the H.E.S.S. 
data on the time scale of only a few minutes, which in itself 
poses a challenging problem. This issue is discussed briefly in 
Section|5] 



1 6b = |Tb ( 1 - /?b cos 6>)] ' , where ySb is the velocity of the emission 
region in c units, T b the bulk Lorentz factor, and 6 the angle between 
the jet axis and the line of sight. 



Throughout the paper, a flat cosmology with Hq = 70 km s 
Mpc~', Qm = 0.3, and = 0.7 is assumed. Unless otherwise 
indicated, all errors are statistical and are given at the lcr confi- 
dence level for one parameter of interest. Dates of observations 
are in general given in the Modified Julian Day (MJD) format. 

2. Observations and data analysis 

2.1. VHE y-ray data 

Due to its high brightness in the VHE range, PKS 2155-304 is 
a regular target for observations with H.E.S.S. and has been fol- 
lowed every year since the start of operations. Scheduled for ob- 
servations towards the end of July 2006, the blazar was found 
in a state of high activity on the first night of the observational 
period, the night of July 25/26 (MJD 53941), with a significance 
of 35cr (standard deviations) for 1.31 hours of live time, which 
is on average 89% of the flux from the Crab nebula as observed 
by H.E.S.S. above 200 GeV On the following night, the source 
was seen with 29.5cr for 1.76 hours of live time. On the night 
of July 27/28 (MJD 53943), the night before the first flare, the 
source was at a lower, but still relatively high flux level, at 12.8<x 
for 1.33 hours of live time. This exceptional activity triggered 
Target of Opportunity observations on the source with SWIFT 
and RXTE; the monitoring however started only one day after 
the major flare of July 28/29 (MJD 53944), where PKS 2155- 
304 was found at 172. 9cr for 1.33 hours of live time. 

The source remained in a relatively high state during the 
whole campaign, with fluxes on most nights exceeding the qui- 
escent state, derived fro m data taken between 2005 and 2007 
dAbramowski et al~ll2010l) . by an order of magnitude. The two 
flares detected on the nights of the 28th and 30th of July 2006 
reached peak fluxes close to two orders of magnitude above the 
quiescent state. The H.E.S.S. light curves for all the nights of 
the 2006 observations are included in Figure Q] together with 
light curves from the other observed wavelength bands, exclud- 
ing the long-term optical and radio data described in Sections l2~3l 
andl24l 

Details of the observations a nd data analysis o f the 2006 
H.E.S.S. data have been given by Abr amowski et al.l d2010l) . To 
extract light curves from the H.E.S.S. data, the VHE y-ray flux 
was integrated above 200 GeV or above the (zenith-angle depen- 
dent) analysis threshold, if the latter was higher. After extraction 
of the spectral shape for each night, integrated fluxes were then 
determined above 300 GeV. 

The 2006 MWL data set includes seven nights with RXTE 
X-ray pointings and eight nights with SWIFT-XRT X-ray 
pointings, which were performed during H.E.S.S. observa- 
tions. Unfortunately, no MWL coverage exists for Flare 1. The 
H.E.S.S. spectra for the nights with X-ray coverage, extracted 
between 0.3 and 2 TeV based on the hypothesis of a power 
law, are shown in Figure |2l with the corresponding photon in- 
dices listed in Table I A. 1 1 in Appendix [A] The spectrum of the 
quiescent state is included for comparison. The spectra have 
been determined from the measured distributions of signals 
from the source and from an off-source region using a for- 
ward folding maximum-li kelihood method (for more details cf. 
Abramowski et al. (2010)). The forward folding method yields 
the best set of parameters for the power law hypothesis with their 
associated errors and the residuals between the measured and ex- 
pected excess in each energy bin. To allow an easier comparison 
with the other wavelength bands, spectral points have been de- 
termined from the residuals. The x 2 values given in Table lATl 
not provided directly by the maximum-likelihood method, have 
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Fig. 1. Light curves in the VHE, X-ray, UV and optical bands during the MWL campaign. Details on the different data sets are 
given in the text. "Flare 1" and "Flare 2" refer to the two exceptional VHE flares seen before and during the MWL campaign. All 
light curves are binned in intervals of 10 minutes, except for the nightly averaged H.E.S.S. fluxes (first panel) and the averaged 
C ANGAROO fluxes (second panel). The published flux from CANGAROO has a threshold of 660 GeV. It has been extrapolated 
down to 300 GeV for comparison with H.E.S.S. The H.E.S.S. points marked in red correspond to nights where simultaneous data 
are available either from RXTE or from SWIFT XRT T he corresponding sim ultaneous points are also marked in red in the light 
curves from the latter. The SWIFT UVOT flux densities dFoschini et al. Il2008l) in the UVW2 and v bands have been multiplied by 
1928 A and 5500 A, respectively. The UVOT error bars include systematic errors. 
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Fig. 2. Left panel: spectra obtained with SWIFT XRT and RXTE during nights of the 2006 MWL campaign. Only nights where 
spectra for both instruments could be extracted are included. For MJD 53946 only the RXTE spectrum from the first pointing and 
for MJD 53947 only the spectrum from the second RXTE pointing have been included to avoid cluttering the plot. The spectra have 
been corrected for Galactic absorption. Right panel: nightly averaged H.E.S.S. spectra for the same selection of nights. Upper limits 
at 99% confidence level are indi cated as arrows. No corre ction for absorption on the extra-galactic background has been applied. 
The quiescent flux is taken from Abramow ski et al.l (l2010t> . The legends in both panels provide the start times of the observations 
in MJD - 53000. 



been determined from power law fits to the spectral points. The 
power law yields a good description for all nights, except for the 
data set fro m Flare 2 (starting at MJ D 53945.87), as discussed 
in depth by lAharonian et aiTI ( 2009bl) . For the night of Flare 2, 
the power law merely provides a rough approximation of the av- 
erage spectrum, which is still useful for a comparison with the 
other nights. 

It should be noted that the H.E.S.S. data during MJD 53949 
were taken with only two telescopes. Nevertheless, the observed 
events were of sufficient quality to be used for the extraction 
of a spectrum for this night, although systematic uncertainties 
are larger than for the other nights. Given the higher threshold 
for observations with only two telescopes (~ 420 GeV in this 
case), the integral flux above 300 GeV for this night had to be 
determined using an extrapolation of the spectrum. Data from 
the three nights MJD 53951 to MJD 53953, where the source 
was in a relatively low state and exposure times were short, have 
been averaged to derive a spectrum with good statistics. 

The nightly averaged VHE flux varies by about an order 
of magnitude for the nights with simultaneous X-ray data. The 
comparison with the spectrum of the quiescent state shows that 



all nights of the 2006 campaign have a significantly elevated flux 
level. 



2.2. X-ray data 

The exceptional activity detected by H.E.S.S. triggered an RXTE 
observation on the source; the monitoring started only one day 
after Flare 1, on MJD 53945. The PCA (lJahodaetal.l fl996) 
units of the RXTE telescope observed the source with expo- 
sures of typically « 700 s per pointing. One pointing was per- 
formed per night, except for two nights (MJD 53946 and 53947) 
with two pointings each. The STANDARD2 data were extracted 
using the HEASOFT 6.5.1 analysis software package provided 
by NASA/GSFC, and filtered using the RXTE Guest Observer 
Facility (GOF) recommended criteria. The spectra were ex- 
tracted using XSPEC v . 12 . 4 . 0, with a fixed colu mn density of 
N H = 1 .7 x 10 20 cm" 2 (iDickev & Lockman II 19901) . A power law 
model was used for spectral fitting. No significant improvement 
was noted when using a broken power law. 

A comparison of the RXTE energy spectra is shown in 
Figure [2] with the results of a power law fit given in Table IA.2I 
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in Appendix lAl The corresponding light curves, binned in time 
intervals of 10 minutes, are included in Figure Q] The lowest 
X-ray fluxes are close to the level detected in the 2003 MWL 
campaign, when the source was in a rather low state and the in- 
tegral flux bet ween 2 and 10 keV was measured at 2.66 x 10" 11 
erg crrT 2 s _1 (Aha ronian et al. 1 l2005bh . 

Simultaneous data from the CHANDRA telescope were only 
available for the night of Flare 2. The decreasing part of an X-ray 
flare is visible during that night, occurring simultaneously with 
the VHE flare. The light curve (integrated flux between 0.5 and 
5 keV) for this night, rebinned in 10 min intervals, is included in 
Figure Q] A detaile d description of tho se data and the analysis is 
published bv lAharonian et al. I (l2009bl) . 

Data taken with XRT onboard SWIFT during the H.E.S.S. 
2006 campaign show an X-ray flare during the night of Flare 
2, which then dec reased by a factor of * 5 in a month 
dFoschini et al. 1 120071) . The light curves of SWIFT XRT and 
SWIFT UVOT are included in Figure [T] rebinned in 10 minute 
bins. SWIFT XRT data are available for seven nights, six out 
of which coincide with nights where RXTE pointings were per- 
formed as well. 

The SWIFT XRT data for the nights during the H.E.S.S. cam- 
paign have been analysed in the present work assuming a col- 
umn density Nh = 1.7 X 10 20 cirT 2 , to be directly comparable 
with the CHANDRA and RXTE data. This analysis has been car- 
ried out with the HEASOFT 6 . 5 standard tools. The data were 
binned requiring a minimum of 20 counts per bin and fitted to 
a single or broken power law. Only points in the intervals from 
0.3 to 0.45 keV and from 0.6 to an upper limit of between 4 
and 7 keV (depending on the data set) were included in the fit 
to suppress known syste matic effects at intermediate energies. 
As iFoschini et al. I (120071) have already pointed out, the spectra 
of several nights are best fit by a broken power law. In the new 
analysis, this was the case for all nights except for MJD 53945, 
53952 and 53953, where a broken power law did not yield a bet- 
ter result than a single power law. During those three nights, the 
duration of the pointings was significantly shorter than for the 
other nights and sufficient photon statistics are only available 
up to about 4 keV. The data recorded during the night of MJD 
53947 do not provide sufficient statistics to extract a spectrum. 
The resulting photon indices and de-absorbed fluxes are listed in 
Table IA.3I in Appendix [A] together with the start times and the 
durations of the pointings. 

2.3. Optical and UV data 

PKS 2155-304 was observed with several telescopes in the op- 
tical and UV range during the H.E.S.S. campaign. Even though 
the coverage of the nights of the campaign is not complete, a 
compilation of the different available data sets provides a good 
picture of the behaviour of the source in these bands. The op- 
tical flux from the host ga laxy is estimated to have an appar- 
ent R magnitude of 15.1 (lAharonian et al. I l2005bl). based on 
measurements in the optical and ne ar infrared (lFalomol [l996: 
iKotilainen. Falomo & Scarpa I Il998h . This contribution is not 
significant for a discussion of the flux variation and is therefore 
neglected in the following. 

For five of the nights with H.E.S.S. observations, data in the 
optical range (V band ) are available from the South African 
Bronberg observatory (Mon ard. B. 1 12007b . Observations were 
made with a 35 cm f/8 "Meade" telescope. Calibration of the 
frames was obtained by calculating magnitude shifts relative to 
two bright isolated stars in the field. The error on the average is 
less than 0.015 mag. The first star was used as a reference and 



its very stable light curve is shown together with the light curve 
obtained from PKS 2155-304 in Figure [3] 

In addition to the observations that were simultaneous with 
the H.E.S.S. campaign, light curves from one night in August 
2006 and from three nights in September 2006 were available 
and have been included in this figure. For all light curves, series 
of six successive data points were averaged to yield a mean value 
and error of the mean over a time period of approximately 151 
seconds. This bin size was chosen, after comparison with the 
signal from the reference star, to average out fluctuations in the 
data induced mostly by variability in the atmosphere. 

Optical and UV data from SWIF T UVOT are also ava ilable 
for the H.E.S.S. observational period (IFoschini et al. 12 007). The 
UV flux during the H.E.S.S. campaign is larger by a factor of 
about 1.5 compared to data from April of 2006. 

Light curves from the Watcher telescope, situated in the 
South African Boyden ObservatorjQ, are shown in Figure [4] for 
the R, V and I optical band. Watcher is a robotically controlled 
40 cm f/14.25 classical Cassegrain telescope. The flux from 
PKS 2155-304 has been calibrated against at least two and up 
to five reference stars in the field of view. The uncertainty in 
each point of the light curve is estimated as the RMS of the fluc- 
tuations around the average fluxes from the reference stars. A 
rotation of the CCD camera by 180° close to midnight during 
each night led to a systematic variation in the brightness of the 
reference stars by up to 6%. A residual effect on the light curves 
of PKS 2155-304 could not be completely corrected in the anal- 
ysis, which led to discrepancies with the Bronberg light curves 
for Watcher data taken after the rotation. It was decided to ex- 
clude all data points for each night after the rotation occurred. 
After this correction, the Watcher and Bronberg light curves are 
found to be in very good agreement. 

Coverage of the source over a longer time period was car- 
ried out with the robotically controlled ROTSE-IIIc telescope, 
located on the H.E.S.S. site. The resulting light curve for all data 
taken in 2006 is shown in Figure [5] The complete data set of 
available ROTSE data from April 2004 to November 2008 is in- 
cluded in Figure [8] and will be discussed in Section [3] ROTSE 
has a wide field of view (1.85° x 1.85°) and is operated without 
filters. A relative R magnitude is derived by comparison of the 
in strumental magnitude with the USNO catalogue as described 
bv lAkerlof et al. ( 20001). A de tailed description of the system is 
given bv lAkerlof et al. I (l2003[) and more informat ion on the data 
analysis can be found in lAharonian et alH (l2005bl) . During the 
period where simultaneous Watcher data are available, the form 
of the ROTSE and Watcher light curves are in good agreement, 
although the normalization of the ROTSE data is lower by about 
0.2 mag due to known systematic biases in the ROTSE data. In 
2003, during a low state of the source, an optical flux between 
13.3 and 13.7 in relative R magnitude was measured by ROTSE. 
Figure shows a higher optical activity and a clear increase of 
the optical flux in 2006. 

2.4. Radio data 

Radio data from three observatories are available for the period 
of the 2006 H.E.S.S. campaign, as well as for a longer interval of 
several months and years before and after the campaign. A long- 
term light curve including measurements at different frequencies 
is shown in Figure [6] 

A monitoring programme with the Nangay Radio Telescope 
(NRT) of extragalactic sources observed with H.E.S.S. is in 



website: http://www.assabfn.co.za/friendsofboyden/boyden.htm 
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Fig. 3. From the upper left to the lower right: PKS 2155-304 light curves from Bronberg optical data (black points) for the five nights 
of observation that coincided with the H.E.S.S. campaign in 2006 (MJD 53945, 53946, 53947, 53951, 53952), for the night MJD 
53954 and for three nights in September (MJD 53987 to 53989). The flat light curve of a reference star is shown for comparison 
(light blue points). For the nights with simultaneous H.E.S.S. coverage, the VHE light curve is shown for comparison (red points, 
y-axis on the right). The data have not been corrected for Galactic absorption. 



place since 2001 . The NRT consists of a single-dish antenna with 
a collecting area of 200 x 34.56 m 2 . Observations of PKS 2155- 
304 were made at a wavelength of 11 cm, wi th a half-power 
beam -width (HPBW) of ~ 2' (RA)x 10' (Dec.) dTheureau et al. I 
120071) . and at 21 cm, with HPBW ~ 4' x 22'. Each observation 
was based on between 2 and 36 one-minute drift scans at fre- 
quencies of 2679 and 269 1 MHz, as well as 1407 and 1420 MHz. 
The flux was calibrated against the stable radio source 4C 55.16 
for the period from MJD 53558 to MJD 54100, covering the 
H.E.S.S. campaign. The calibrators 3C 123, 3C 161, 3C 286 and 
3C 295 were used for the data analysis before and after this pe- 
riod. 

The Hartebeesthoek Radio Astronomy Observatory 
(HartRAO), located in South Africa, is equipped with a 26 m 
diameter sing le dishQ Observations with HartRAO started on 
July 28, after an alert from the H.E.S.S. collaboration. Between 
3 and 5 scans were performed each day for the following two 
weeks. Flux measurements at 6 cm (HPBW ~ 10') and 3.5 
cm (HPBW ~ 6') were carried out with dual-beam receivers. 
Continuous flux measurements of PKS 2155-304 used the drift 
scan method, with scans offset at the half -power points North 
and South to check pointing accuracy. The radio galaxy 3C 123 



was used as a calibrator. Its flux density was seen to vary by less 
than 2%. 

One or two daily scans were carried out in the months af- 
ter the campaign up to MJD 54741 (October 1st, 2008). Scans 
where the two receivers differ by more than 10% were rejected. 
The daily fluxes from both receivers were used to calculate aver- 
age flux densities and to estimate the uncertainty. A systematic 
error of 6% - 7% was estimated to account for uncertainties in 
the data acquisition and in the data reduction procedure. 

The Australian Telescope Compact Array (ATCA) is located 
at the Narrabri Observatory and consists of six 22 meter anten- 
nasS ATCA data at 3 cm (HPBW ~ 5'), 6 cm (HPBW ~ 9'), 
13 cm (HPBW ~ 20') and 20 cm (HPBW ~ 31') are included 
in Figure [6] Data were taken in the snapshot observing mode, 
which does not provide angular resolution. Only averages of 
observations from July 28 to August 4 are available. The ob- 
servations start at about 14:00 UT and do not overlap with the 
H.E.S.S. observational window. No significant nightly flux vari- 
ation was found, but the overall flux density is higher than that 
detected in previous observations from mid 1997 to mid 2000. 



website: http://www.hartrao.ac.za 



website : http ://w w w. narrabri . atnf . csiro.au 
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Fig. 4. Light curves of Watcher optical data in the I band (upper 
points; in black), R band (points in the middle; in red) and V 
band (lower points; in green). The arrow indicates the time in- 
terval of the MWL campaign. The data have not been corrected 
for Galactic absorption. 
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Fig. 5. PKS 2155-304 light curve observed with the ROTSE 
telescope in 2006 (R band magnitude). The arrows indicate the 
time intervals of the MWL campaign and of coverage with the 
Watcher telescope. The data have not been corrected for Galactic 
absorption. 



3. MWL spectral and flux evolution 

3.1. Flux evolution in different energy bands and correlations 

The VHE and X-ray light curves and a sample of the UV and op- 
tical light curves available from MJD 53945 to MJD 53953 are 
compiled in Figure Q] All data have been binned in 10 minute 
bins, with the exception of the light curve in the uppermost 
panel, which provides the nightly averaged flux observed with 
the H.E.S.S. telescopes and the average fluxes detected with the 
CANGAROO experi ment, which correspon d to a livetime of 
typically a few hours (Sakamoto et al. 2008). Only the nightly 
averaged three-fold triggers are shown for CANGAROO. To 



allow a more direct comparison, the published fluxes from 
CANGAROO, measured above 660 GeV, have been extrapolated 
down to the same threshold as the H.E.S.S. integrated fluxes (300 
GeV). The photon indices determined with H.E.S.S. from the 
data sets closest in time to the CANGAROO measurements were 
used in this extrapolation. The CANGAROO data complement 
the H.E.S.S. data and indicate that the source was in a low state 
just before Flare 2 occurred. 

A first look at the other wavelength bands reveals a high flux 
in the X-ray band during the night of Flare 2. Observations with 
CHANDRA started only after the peak in the VHE flux and show 
a decreasing flare simultaneous to Flare 2. Data from RXTE and 
SWIFT XRT, taken over much shorter time intervals than the 
data from Chandra, reveal an X-ray flux during this night that 
is higher than on the previous and on the following night and 
that is on the decrease. The optical flux is increasing during that 
night (MJD 53946), but is lower than during several of the other 
nights. 

For the SWIFT XRT and RXTE data sets, light curves have 
been extracted between 0.3 and 4 keV and between 2 and 10 keV, 
respectively, in intervals of 10 minutes. For statistical reasons, 
nightly averaged spectra were used to calculate the absolute val- 
ues of the de-absorbed fluxes. It is assumed that intra-night vari- 
ations in the photon index are negligible for the determintation 
of the integrated fluxes. The XRT and RXTE light curves can be 
seen to follow the same evolution. The X-ray flux is in its low- 
est state on the first night (MJD 53945) and shows its largest 
increase on the night of Flare 2, although X-ray data are only 
available after the peak observed in the VHE range. In the nights 
following Flare 2, the X-ray flux is decreasing more slowly than 
the VHE flux, while staying always above the low level of the 
first night. 

In the UV and optical bands, the flux densities have been 
multiplied with the effective frequencies of the corresponding 
bands to allow an easier comparison of the energy output at dif- 
ferent wavelengths. (It should be noted that the optical filter in 
the V-band used by SWIFT UVOT is different from the one used 
with the Bronberg and Watcher telescopes.) The fluxes included 
in Figure[T]have not been de-reddened. The UV and optical flux 
increases during and after the night of Flare 2 and stays above 
its initial level in the following two or three nights, but the small 
amplitude of this increase is comparable to variations observed 
on other nights. 

Peak to valley flux variations of about a factor of four can 
be seen in the X-ray range over the selected period, whereas 
the VHE flux is increasing by up to two orders of magnitude. 
Variations in the optical range are much smaller, of the order of 
30% in the V band. 



3.1 .1 . VHE and X-ray flux evolution 

It has previously been shown that the VHE y-ray flux corre- 
lates with almost the third po wer of the X-ray flux o bserved by 
CHANDRA during Flare 2 (lAharonian et al. ll2009bl) . When ex- 
amining the evolution of the X-ray and VHE flux (Fig.|7]i during 
the other nights of the campaign, one can see that, also on a 
longer time scale, variations in the VHE flux correlate to varia- 
tions in the X-ray flux observed with RXTE and SWIFT XRT. 

When evaluating the correlation between RXTE and H.E.S.S. 
over the whole campaign (using a 10 to 90 minute bin- 
ning, see Fig. |7), a correlation coefficient of 0.75 was 



5 Correlation coefficients here and in the following sections corre- 
spond to the "Pearson r". 
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Fig. 6. Long term radio light curve with data from the NRT, ATCA and HartRAO, from July 2004 to November 2008. The lower 
panel is a zoom on the period around the H.E.S.S. campaign. HartRAO error bars include an estimate of the systematic errors. The 
double arrow marked "MWL" in the lower panel indicates the time period of the MWL campaign. 



found, c orresponding to a cha nce probability of less than 
0.001% dPugh & Winslow 111966b . The comparison of H.E.S.S. 
and XRT data yields a correlation coefficient of 0.72 with a 
chance probability of less than 0.03%. When ignoring all data 
points from the night of Flare 2, the remaining data still yield 
correlation coefficients of 0.61 (RXTE) and 0.54 (XRT) with 
chance probabilities of less than 0.5% and less than 5%, respec- 
tively. 

Power law fits to subsets of simultaneous data points are 
shown in Fig. [7] They do not provide a good description of 
the data in statistical terms, hinting at a more complicated be- 
haviour, but they show a general trend of the correlation pattern. 
For both data sets, data points from nights with relatively high 
fluxes have been fitted to illustrate this trend. Arbitrary values 
of log F(2-10keV) > -10.2 log(erg cirT2 s"l) and log F(0.3-4 
keV) > -9.65 log(erg cirT2 s~l) have been chosen to define the 
subset of data with relatively high fluxes. The fits to the "high 
flux" points have very steep slopes, with exponents between 4 
and 5 (dashed lines). When removing all points from the night 
of Flare 2, the data from the "high flux" states can be fitted with 
exponents between 2.5 and 3.5 (solid lines) , still presenting a 
steep correlation between the VHE and X-ray bands. 

The power law fits indicate that the observed correlation be- 
tween VHE and X-ray fluxes varies over the period of the MWL 
campaign. At high flux levels, the VHE flux varies as a high 
power of the X-ray flux variation. The correlation becomes less 
steep with decreasing flux and there is no clear correlation for 
the nights with the lowest flux levels. This indicates a different 
behaviour of the emission depending on the flux level. It is also 
coherent with a recent joint H.E . S.S./F ermi//?XrE'/ATOM cam- 
paign in 2008 ( Aha ronian et al. 2009a), where the source was 



found in a relatively low state and no correlation between the X- 
ray and VHE band was detected, similar to the 2003 low state. 
Taken together with the behaviour described here, this might in- 
dicate that a correlation between these bands occurs only if the 
source is in a high state. 

It should be noted that evidence of a difference in the 
behaviour of the emission spectrum of t his source between 
low- a nd high-flux states was also found by lAbramowski et al. I 
(2010) concerning the evolution of the VHE spectral index with 
the VHE flux level. 

A correlation between the X-ray and VHE fluxes is nat- 
urally accounted for in the SSC framework and has been 
detected for example during high-flux state s of Mrk 421 
dKrawczvnski et aTTT 2001 : Fossati et al. 2004: iTanihata et aF 
1 2004 1 : iBIazeiowski et al 1 120051) and Mrk 501 (ICatanese et al 
1997] iPian et al. 1119981; iDiannati-Atai et al. Ill999t lAlbert et al. 
2007). An extension to more than one emission zone has been 
proposed to account for the observed steep co rrelation between 
the X- ray and VHE y-ray flux during Flare 2 (Ahar onian et al. I 
2009b). An application of multi-zone SSC models to Flare 2 will 
be demonstrated in Section [4] Such a scenario could be appli- 
cable more generally to explain the changing behaviour of the 
source between high- and low-flux states, as will be discussed in 
Section|5] 



3.1.2. Optical flux evolution 

As can be seen from the ROTSE data in Fig. [5] the optical 
flux from the source was at a relatively high level during the 
2006 MWL cam paign, compared to the l ow state observed with 
ROTSE in 2003 d Aharonian et al. Il2005bl) . This is also visible in 
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Fig. 7. Correlation between the integrated VHE flux and the integrated RXTE (left) and SWIFT XRT fluxes (right). The labels 
indicate the dates of observations (MJD - 53000). The points represent averages over simultaneous time intervals of 10 minutes 
(except for the H.E.S.S./RXTE point from MJD 53949 and the H.E.S.S./XRT point from MJD 53952 that correspond to a 30 
minute time window, as well as the H.E.S.S./XRT point from MJD 53953, corresponding to a 90 minute interval). The statistical 
uncertainties on the logarithmic fluxes are approximated with symmetric error bars. The two lines in each panel show power law 
fits to subsets of the data: the dashed lines include all "high flux" data, while the fits marked with solid lines include all "high flux" 
data except for the data from the night of Flare 2 (MJD 53946). 



the ROTSE long-term light curve (Fig. [HJ and is confirmed by a 
compilation of archival data and by obs ervations with the robotic 
60 cm telescope REM in 2005 (iDolcini et al. 112007b , indicating 
typical fluxes from PKS 2155-304 in the V band to lie between 
about 16.5 and 26 mJy. During the 2006 campaign, fluxes of 
roughly 27 to 33 mJy were observed by the Bronberg observa- 
tory and SWIFT UVOT 

On intra-night time scales, no exceptional variations are ob- 
served in the optical flux during the active VHE state of the 
source. In the Bronberg, Watcher and ROTSE data, intra-night 
variations of «0.1 mag can be seen for several nights during the 
MWL campaign and also during the following months. The min- 
imum time scale for these variations is of the order of 1 hour, 
corresponding to the width of the small peak in the Bronberg 
light curve around MJD 53987.85. The amplitude of the ob- 
served variability is similar to observations during the 2003 low 
state. 

During the night of Flare 2, the Bronberg data show an in- 
crease of »0.1 mag (see Fig. [3}. Over the following nights, the 
optical and UV flux rise to a local maximum that occurs proba- 
bly two or three nights after the flare, judging from a comparison 
of the available Bronberg, Watcher and SWIFT data (see Fig. [T}. 
Several similar episodes of flux increases over a few nights are 
seen in the Watcher and ROTSE data and seem to b e a general 
featur e of the optical flux for this source (see e.g. iRvle et al~l 
2006). 



The correlation coefficients between the integrated H.E.S.S. 
and Bronberg fluxes for the different nights, taken at face 
value, indicate a positive correlation for MJD 53945 (0.76 for 
34 points) and for MJD 53952 (0.98 for 6 points), an anti- 
correlation for MJD 53946 (-0.74 for 41 points), no correlation 
for the nights of MJD 53947 and 53951. The ubiquitous varia- 
tions in the optical band lead to random instants of correlations 
and anti-correlations with the VHE band. When taking into ac- 
count the overall behaviour of the optical flux for all the nights 
where detailed optical light curves were available, there is no 
compelling evidence of a direct correlation of the optical band 
with the VHE band. 

However, on a time scale of several weeks, the average 
optical flux was rising, as can be seen in the Watcher and 
ROTSE data and by comparing the three September nights in 
the Bronberg data to the observations from July and August. 
During the Watcher observation period, the flux increased by 
»0.6 mag, with the fastest increase being ^QA mag in less than 
10 days, occurring after the MW L campaign, when the VH E flux 
was seen in a rather low state dAbramowski et alTl 12010). The 
ROTSE data, which cover a longer time span, show a fluctuation 
of * 1.2 mag between the lowest and highest fluxes in 2006. The 
source was observed with REM from August 23 (MJD 53970) 
on and an increase to a particularly bright state on October 17 
(MJD 54025) of V = 12 mag was found, followed by a rapid de- 
crease. Variations in the V band of abou t 0.6 mag were detected 
by SWIFT UVOT dFoschini et al. 120081) . The long term rise seen 
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by ROTSE, Watcher, REM and Bronberg is not reflected in the 
VHE y-ray data, but suggests a relation between high optical 
states and flaring VHE activity. 

During the 2008 MWL campaign on PKS 2155-304, some 
evidence was seen of a direct correlation between the VHE and 
optical fluxes on a time scale of several nights. If such a correla- 
tion during the low state of the source is confirmed by future ob- 
servations, it would be another piece of evidence of a difference 
in the emission characteristics between the high and low states 
of the source. Stronger statements on this issue require long-term 
MWL monitoring of the source, not only covering flaring states, 
but also regularly sampling low flux states. 

3.1.3. Radio flux evolution 

The evolution of the available radio data (Hig. [6j shows that the 
source was already in an active state in the radio band during 
the H.E.S.S. campaign and exhibited a significant increase in the 
flux level about one and a half months after the VHE flares oc- 
curred. No significant variation was seen between the nightly 
ATCA fluxes during the MWL campaign, which is confirmed by 
the constant flux in the HartRAO data. The average flux mea- 
sured by the NRT during the campaign (0.49+0.04 Jy at 1 1 cm) 
is significantly higher than the flux measured during the 2003 
campaign (0.30+0.01 Jy at 11 cm). This is consistent with the 
ATCA flux being higher than previous measurements found in 
the ATCA database. 

In Figure [8] the long-term light curve in the radio band 
(NRT, HartRAO) is compared with the optical light curve from 
ROTSE and with data from H.E.S.S., taken from 2005 to 2007. 
The rise in the radio flux in the months around MJD 54000 
(September 2006) by roughly 50%, observed by both the NRT 
and HartRAO telescopes, coincides with a comparable increase 
in the optical flux during that time period. It can also be seen that 
both the radio and optical fluxes were at a lower level around 
MJD 53200 (July 2004) and decline again around MJD 54400 
(October 2007), even though in this case the relative change in 
the amplitude is not the same. Both bands exhibit a relatively 
high flux level that increases still during or shortly after the 
VHE high state. Although the optical and radio flux evolution 
shows some differences, especially around MJD 53300 and MJD 
54400, which might indicate a delay between the two bands that 
is difficult to quantify due to the sparse data coverage, an overall 
correlation is clearly present. Over the whole period presented 
here, a correlation coefficient of 0.29 for 128 points is derived, 
corresponding to a chance probability of 0.1%. 

More recent data from the Effelsberg radio telescope show 
that the radio flux from PKS 215 5-304 was continuous ly de- 
creasing further in 2009 and 2010 ( Angelaki s et al. I l2010h . 

3.2. Spectral variability 

Spectral hardening with an increasing flux level is usually ex- 
plained with the injection of highly energetic particles into the 
emission zone or with rapid particle acceleration. The observa- 
tion of spectral hardening at both X-ray and TeV en ergies, ob- 
served for example during th e 1997 flare of Mrk 501 dPian et al. I 
1998l lTavecchioet al. 12 001). is evidence of a common origin of 
the variability in those bands, as expected in the SSC framework. 

In the case of PKS 2155-304, spectral hardening with in- 
creasing flux h ad been found in the RX TE data during the low 
state in 2003 dAharonian et al. I l2005bh and i n the RXTE and 
SWIFT XRT data during the low state in 2008 dAharonian et al. I 



2009a]), but no significant spectral variability was seen in the 
VHE band at that time. On the other hand, the H.E.S.S. data 
taken during the four nights with the highest fluxes of the 2006 
high state (MJD 53944 to MJD 53947) show a clear hardening 
of the VHE spectrum with flux increase. Finally, the analysis of 
H.E.S.S. data from 2005 to 2007 confirms that the source does 
not ex hibit the same spectral b ehaviour at high- and low-flux 
states dAbramowski et al7ll2010l) . 

The MWL spectral evolution during the 2006 campaign has 
so far only been studied on the short intra-night time scales of 
Flare 2, wh ere spectral hardening is seen in both the X-ray and 
VHE band dAharonian et al. |[2009bl) . The spectral behaviour of 
nightly averaged VHE fluxes over the whole 2006 MWL cam- 
paign, in comparison with X-ray fluxes from the same period, is 
shown in Fig. [9] The RXTE data show no indication for a corre- 
lation between the photon index and the hard X-ray flux (inte- 
grated between 2 and 10 keV). Spectral hardening can be seen 
in the SWIFT XRT data in the soft photon flux below the break 
energy of «1 keV (cf. Tab. IA.3I in Section l2~2l ). Even though 
the overall change in the photon index is very modest (from 2.84 
to 2.37), the correlation coefficient of 0.93 indicates a very sig- 
nificant effect. This correlation remains significant, albeit with a 
smaller coefficient of 0.86, if one imposes an ad hoc break in the 
spectrum at 1 keV for the data from MJD 53945 and MJD 53952 
to be able to apply a broken power law for most of the nights (cf. 
Tab. IA.3b . Above the break energy, no significant variation in the 
photon index is seen in the XRT data. Thus the 2006 high state of 
the source seems to be characterized by spectral hardening only 
in the soft X-ray range. Faster cooling at higher energies might 
suppress the effect of spectral hardening in the hard X-ray band. 

During the nights of the MWL campaign, the spectral in- 
dex of the nightly averaged H.E.S.S. spectra changed by »0.8. 
For this selection of nights, the overall correlation coefficient 
of 0.59 is not significant, but the spectral evolution gives a 
hint of a pattern more complex than a linear correlation. The 
H.E.S.S. data show some indication for a counter-clockwise loop 
in the last four data points. On a shorter time-scale, a similar 
behaviour has already been seen in the VHE ra nge with the 
MAG IC telescope during a flare from Mrk 50l dAlbert et al. I 
120071) . Features in the form of a hysteresis are generally ex- 
pected in X-ray flares that are su bject to synchrotron cooling 
dKirk. Rieg er & Mastichiadis 1 1998b and have been observed for 
example with SWIFT X RT during flares in Mrk 421 in 2006 
dTramacere et al. 1120091) . or in PK S 2155-304, with th e LAC on 
board the Ginga satellite in 1993 dSembav et al. II 19931) and with 
ASCA in 2000 dKataoka et a l. 2000). Evidence of loops have 
also been found recently in the high energy range (above 100 
MeV) with the Ferm i telescope in the case of the flat spectrum 
radio quasar 3C 273 (lAbdo et al. II2010I) . 

In the X-ray range, the emergence of clockwise and 
counter-clockwise loops is well understood. According to 
iKirk. Rieger & Mastichiadisl dl998l) . a clockwise loop can arise 
when cooling is important, but the acceleration time scale is 
much faster than the cooling time scale, which is the case if 
the system is observed below the synchrotron peak. A counter- 
clockwise loop indicates that the acceleration and cooling rates 
are comparable, indicating that the system is observed closer to 
the maximum frequency. In this case, the gradual acceleration 
of particles leads to increased emission first at lower and then 
at higher energies. In the SWIFT XRT data set presented here, 
no conclusion can be drawn on the existence of such features, 
due to the low statistics and the missing spectral information for 
MJD 53947. 
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Fig. 8. Long-term light curves from VHE y-ray measurements (upper panel: H.E.S.S.), optical data (middle panel: ROTSE) and 
radio data (lower panel: NRT and HartRAO). Please note that the y-axis is truncated for the long-term H.E.S.S. light curve due to 
the large dynamic range caused by the flares. A zoom on the period covering the two VHE flares is shown in an inset in the upper 
panel. 



In the SSC framework, where the (VHE) y-ray emission is 
strongly connected with the synchrotron emission, loops in the 
VHE data could reflect related features in the X-ray range. If 
such patterns, as suggested by the present data, were confirmed 
in the VHE range by future instruments with high er sensitiv- 
ity, su ch as the Cherenkov Telescope Array dCTA Consortium! 
2010), this could provide important additional constraints on the 
validity of the SSC scenario. 

4. Modelling of the MWL data on different time 
scales 

The MWL coverage of the source over several nights makes it 
possible to investigate the spectral energy distribution (SED) and 
its evolution on two different time scales. Night- to-night vari- 
ations can be evaluated with stationary models for non-flaring 
nights. Conversely, the night of Flare 2 allows the study of the 
evolution of the SED on intra-night time scales. Interpreting the 
variability observed during this night requires time dependent 
modelling. 

In Section |4~T1 a conventional stationary SSC model is used 
to describe the nightly averaged SEDs, while the flux evolu- 
tion during Flare 2 is illustrated using different examples for 
time-dependent SSC scenarios in Section 14.21 Basic consider- 



ations on the requi rements of a time-depen dent scenario for the 
night of Flare 2 ( Ahar onian et al. 1 1200 9b) challenge the com- 
mon scenarios for VHE blazar emission and show that standard 
one-zone SSC models cannot account for the detected SED evo- 
lution. Here, three different approaches to time-dependent SSC 
modelling with increasing complexity are discussed, which are 
representative of the most common options available in the liter- 
ature beyond one-zone models. The first scenario is a basic de- 
scription of an extended VHE-emitting source, well adapted to 
fast variability processes, focusing on high-frequency emission. 
The second one considers a two-component jet to reproduce the 
SED and its evolution. The third one proposes a stratified jet with 
time-dependent particle injection and acceleration. The first two 
scenarios assume that the initial particle acceleration is due to 
first order Fermi processes, while the third one favors a second 
order Fermi mechanism. 

The intention of the time-dependent modelling of Flare 2 
presented in this section is to illustrate how different SSC ap- 
proaches can arrive at a good characterization of the dominant 
components in the flux variation and what constraints the data 
provide in general on physical source parameters. These mod- 
els present three possible ways of describing the data set and do 
by no means exclude alternative descriptions. A description of 
the very rapid variability observed during Flare 1 and Flare 2 
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Fig. 9. Variation of the photon index with the integrated flux. The 
upper panel shows the variation in the soft band of the SWIFT 
XRT data, the middle panel in the RXTE data, and the lower 
panel in the H.E.S.S. data. Labels indicate the date in MJD- 
53000. The XRT points marked with filled circles show the pho- 
ton index below the break for nights where a broken power law 
yields the best result and the photon index for a simple power 
law for the remaining nights. For MJD 53945 and 53952, solu- 
tions for the XRT adopting an ad hoc break at 1 keV are shown 
with open circles. For more details, see Tab. IA.3I The H.E.S.S. 
point with the label "951" refers to an average over data taken 
from MJD 53951 to 53953. 



dAbramowski et al71l2010h . on time scales down to a few min- 
utes, is beyond the scope of the models presented here, which 
focus on the general trend of the flux evolution during the night 
of Flare 2. 

For all the models discussed in this section, absorption of 
the VHE y-rays b y the Extragalactic Backg round Light (EBL) 
is treated following Aharonian et al. (2009a). 



4.1. Stationary modelling : one-zone SSC model 

The average SEDs for the different nights of the MWL cam- 
paign, except for the night of Flare 2, can be described with a 
stationary homogeneous one-zone SSC model. The X-ray and 
VHE bumps are interpreted respectively as the synchrotron and 
IC emission from a single population of relativistic electrons in a 
spherical plasma blob inside the jet. The host galaxy and the ex- 
tended jet contribute to the SED at lower frequencies. The blob 
is characterized by it s radius r^, Doppler factor 6b a nd tangled 
magnetic field B (see lKatarzvnski. Sol & Kusll200lh . The elec- 
tron energy distribution is assumed to follow a broken power 
law with indices n\, «2, minimum and maximum energies cor- 



responding to Lorentz factors y m \ n and y max and a break energy 
corresponding to y^. The parameter K\ determines the normal- 
ization of the broken power law at a Lorentz factor of 1 . This 
type of SSC model has typically eight main free parameters and 
thus requires high quality data over a large spectral range to be 
well constrained. 

The available data on the nights of the MWL campaign do 
not provide enough constraints on the SEDs to narrow down the 
parameter space even of such a stationary model. Therefore, a 
set of paramet ers derived from t he 2003 H.E.S.S. campaign (cf. 
"Model 2" in Aharonia n et al. ] (|2005b)) was used as a starting 
point for the modelling of the 2006 data of the source, in or- 
der to compare the results for 2003 and 2006. These parameters 
have been adapted to the nightly averages of the 2006 campaign 
by applying only minimal changes to the initial values. A set 
of parameters that reproduces the average SED of the relatively 
low state on MJD 53945 is presented in Table Q] Even though 
the VHE y-ray flux during the whole 2006 observation period 
is higher, the SED of MJD 53945 (Fig.[l0b is not very different 
from the one recorded in 2003. For the same Doppler factor and 
magnetic field strength, a slight increase of the size and density 
of the emitting blob and the spectral index «2 provides a good 
match of the SED observed during MJD 53945. 
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Table 1. Parameters used for stationary SSC modelling of the 
nightly average spectrum of MJD 53945. A description of the 
parameters is given in the text. 



The difference between n \ and «2 is larger than what would 
be expected for a pure synchrotron cooling break in the electron 
spectrum; the latter does not yield a satisfactory presentation of 
the data. The broken power law here is a simple parameterization 
of the electron spectrum, which may be interpreted as a sum of 
several effects that are not explicitly modelled, e.g. synchrotron 
and IC cooling, adiabatic expansion, energy-dependent acceler- 
ation and escape times, inhomogeneous particle distribution, etc. 

Given the absence of direct correlations between the high- 
energy (VHE and X-ray) and low-energy (optical and radio) 
bands in the data, the radio component, and the main part of the 
optical emission is ascribed to a different origin than the high- 
energy bumps. It can be modelled as the emission of the ex- 
tended jet (see lAharonian et al. Il2005bn . with an additional con- 
tamination by the host galaxy. A relatively low electron density 
and magnetic field in the extended jet results in a synchrotron 
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Fig. 10. SED for a night with relatively low VHE and X-ray 
flux (MJD 53945) with a stationary SSC model described in 
the text. Data taken during this night are shown from H.E.S.S. 
(red squares), RXTE (magenta circles), SWIFTXRT (green trian- 
gles), SWIFT UVOT (red triangles) and Bronberg (black circle). 
The black squares are the average ATCA points from the period 
of the MWL campaign. The model is presented as a solid line, 
giving the overall SED, including absorption by the EBL. The 
dashed line in the VHE range is the de-absorbed spectrum. The 
H.E.S.S. data points represent the measured (i.e. absorbed) flux. 
The extended jet emission (dotted line) and host galaxy contri- 
bution (brown data points) are also shown at low frequencies. 
Dash-dotted and dashed lines indicate the synchrotron and SSC 
contributions from the blob, respectively. 



emission peak from the radio to the UV range that is accompa- 
nied only by a negligible IC component at X-ray energies. For 
the extended jet shown in Fig. [10] a jet opening angle and angle 
to the line of sight of 2 degrees have been assumed. The other 
main parameters, such as the radius rj et , length Zj et , the initial 
magnetic field Bj et and electron distribution, are listed in TableQ] 
A more detailed description of all the parameters in this inho- 
mogen eous jet model can be found in Kat arzvnski. Sol & Kus I 
d200lh . 

The contribution of the host galaxy, as shown in 
Fig. [TUl has been deduced from the magnitudes given in 
Kotilaine n. Falomo & Scarpa] (Il998l) . assuming a low redshift 
solar metallicity elliptical galaxy of age equal to 13 Gyrs 
R-H=2.4), corresponding t o a m ass of 5x10" solar masses 
Fioc & Rocca-Volmerangd 1 19971) . The galaxy is clearly not 
contributing much in the optical, even at a state of low activity. 

One can find good solutions for all the remaining nights 
(MJD 53947 to 53950) after Flare 2, with the same bulk Doppler 
factor and magnetic field and a similar electron distribution, 
while only changing the radius of the blob (by less than « 40%), 
the normalization of the electron spectrum (by up to a factor of 
2.4) and «2 (change in index of up to 0.2). The indicator for 
the electron number N tot can be kept constant during these ad- 
justments. Such variations could arise if the shock front passed 
regions of different size and electron density, or for contraction 
and expansion of the emitting blob. 

It should be stressed that the solutions found here are not 
unique, due to the insufficient information from the MWL ob- 
servations. Recent data from the Fermi satellite have shown the 
importance of covering the high-energy y-ray range below VHE 
energies in this respect. 



4.2. Time-dependent modelling : Two-zone and stratified jet 
SSC models 

4.2.1. Model 1 : Light crossing time effects 

A relatively simple app r oach is given by a model based on 
Chiaberge & Ghise llini I d!999l) . A short description can be 
found in Appendix IB. II This model has already been success- 
fully applied to Flare 1 dKatarzvnski et al. II2008I) . but the X-ray 
light curve measured by CHANDRA during Flare 2 adds con- 
straints which were not available for the analysis of Flare 1 . The 
realization of Model 1 that is in best agreement with the data of 
MJD 53946 is shown in Figure[TT]and a list of the corresponding 
parameters is given in Table lB.fl in Appendix lB.il 

Two independent sources were necessary to explain the ac- 
tivity at high energies observed during that night. The first source 
has an extension of 4 x 10 16 cm and an electron density of 
1.6 x 10 4 crrT 3 and leads to a long-term component with a time 
scale below two days that is visible in the X-ray flux and is con- 
sistent with the observations by RXTE and SWIFT XRT on the 
previous and following nights. This source dominates the syn- 
chrotron emission due to its large size, but does not contribute 
significantly to the IC emission due to its relatively low electron 
density. The second source is much smaller (2.35 x 10 15 cm) and 
denser by a factor of 100, with a larger 6b than for the first source. 
This component dominates the IC emission, while contributing 
little to the synchrotron flux. It is responsible for the rapid vari- 
ability of the order of an hour that characterizes the main VHE 
flare. The contributions of the two components to the integrated 
X-ray and VHE flux are shown in the middle and lower panel of 
Figure [TT1 

In this model, the adjustment of few main parameters — 
namely the extension of the emission region, magnetic field and 
density of the injected electrons and the bulk Doppler factor of 
each source — provides a satisfactory description of the main 
components of the observed high-energy data during the two 
VHE flares. The parameters found for Flare 2 are comparable 
to those for Flare 1, with identical <5t,. Additional small and 
dense sources can be added to account for the structure visi- 
ble on smaller time scale s, as has been shown for Flare 1 by 
Katarzvns ki et al. I (120081) . The long term component found in 
Model 1 does not extend over the whole duration of the MWL 
campaign; it is constrained by X-ray data from the night before 
and after the night of Flare 2. Additional large emission zones 
would need to be added to reproduce the behaviour of the X-ray 
flux detected by RXTE and SWIFT over several nights. The op- 
tical and radio emission is ascribed to an additional component, 
e.g. the extended jet, which is not included in this model. Thus 
the low frequency points only serve as upper limits for the SED 
derived with Model 1 . 

4.2.2. Model 2 : Blob-in-jet scenario 

A mor e refined approach proposed by Kata rzvnski. Sol & Kus I 
(2003) is aimed at providing a complete description of the SEDs 
of blazars and of their variability from radio to VHE energies. A 
short description can be found in Appendix lB.2l A first applica- 
tion o f this mod el to Flare 2 ha s been presented bv lLenain et al. I 
(120081) (see also lLenain I d2009ln . 

The application of the blob-in-jet model to Flare 2 is shown 
in Figure [12] with the corresponding parameters in Table IB .2 1 
in Appendix IB. 21 The SED shows the spectra derived from 
three different H.E.S.S. data acquisition periods during that 
night, each of a duration of about 28 minutes, together with 
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Fig. 11. Time dependent SSC model (Model 1) for the night of MJD 53946. The upper panel shows the modelled SED at the 
maximum flux level and at a low level. Snapshots of the CHANDRA spectrum and the H.E.S.S. spectrum are shown for a high state 
(red squares) and for a low state (blue squares). The two flux states that were chosen for the modelled SEDs correspond to the red and 
blue circles in the middle and bottom panel, where the H.E.S.S. and CHANDRA light curves are shown together with the predicted 
time evolution of the model. The thin line in the middle and bottom panel indicates the long-term component that dominates the 
X- ray flux. The two s pectra f rom CHANDRA correspond to the beginning and the plateau of the CHANDRA lightcurve, as described 
bv I Aharonian et al~l d2009bh : the H.E.S.S. spectra correspond to the data acquisition periods closest in time to the modelled flux 
states. 



the area corresponding to the varying spectrum observed with 
CHANDRA. The model describes well the X-ray and VHE spec- 
tra at different flux levels during the night. It also reproduces 
well the main flux evolution at X-ray and VHE energies, as can 
be seen from the comparison between model and observed light 
curves. The observation with CHANDRA began later than the 
H.E.S.S. observation that night, and the rising part of the flare 
is not covered in the X-ray data, so the highest flux state is not 
constrained in the X-ray band. 

Contrary to Model 1, the two high-energy emission regions 
(jet and blob) of Model 2 are not fully independent. Here it is 
assumed that they have the same bulk Doppler factor of 50. 
The blob is denser than the jet by three orders of magnitude, 
with higher magnetic field, higher maximum particle energy and 
smaller initial radial extension. The ratio of the energy density 
in the magnetic field and in the radiation field in the initial emis- 
sion region is very large for the jet (10 3 ), which is the princi- 
pal synchrotron source, and close to unity for the blob, which 
is mainly responsible for the SSC emission. However, the initial 



extension of the blob, its initial electron density and magnetic 
field are close to the values of the small and dense "source 2" in 
Model 1 . The main differences between the two approaches lie 
in the global AGN scenario and in the resulting values for the 
bulk Doppler factors, which are larger in Model 2. This stems 
mainly from the necessity to compensate for adiabatic losses in 
Model 2, which are not taken into account in the simpler Model 
1 . Additional plasma blobs may be introduced to account for a 
finer temporal structure of the light curves. 

The radio and optical emission is reproduced by the addi- 
tion of a large-scale jet with the same parameters as in the sta- 
tionary model (cf. Table Q]), except for a small decrease in the 
magnetic field to B =0.04 G. The optical emission is dominated 
by emission from the large-scale jet, but has also contributions 
from the plasma blob and inner jet. The rapid variation in the 
emission from the blob, which explains the X-ray and VHE flare, 
leads only to a very slight variation in the optical and UV range, 
smaller than the size of the statistical error bars for the nightly 
averaged flux. 
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Fig. 12. Time dependent SSC model (Model 2) for the night of MJD 53946. Left panel: H.E.S.S. spectra at three different states of 
the source are shown. The varying X-ray spectrum seen by the CHANDRA telescope is indicated by a dark grey band. Snapshots of 
the SED resulting from the model for the three different states are presented by solid lines. Optical and UV data from SWIFT UVOT 
and radio data from ATCA are shown as green points. Right panel: The light curves in the VHE band (H.E.S.S., upper panel) and 
in the X-ray band (CHANDRA, lower panel) are shown in units of 1CT 10 erg cm-2 s _1 , together with the predicted time evolution 
of the model. The times corresponding to the three H.E.S.S. spectra in the left panel are marked here by three vertical dash-dotted 
lines. 






In the optical data, variability is seen on different time scales 
and might be due to more than one emission region. The vari- 
ability time scales can be used to constrain the size of the po- 
tential emission regions as R < c ? var 6/(1 + z). For the intra- 
night variability seen in the Bronberg data on the hour scale, this 
would lead to a ratio of source extension and Doppler factor of 
roughly R 5 < 10 14 cm. The derived 5b\ob and R\,\ b have a ratio 
of the same order of magnitude, so that the plasma blob might 
in this model be identified as the emission zone responsible for 
the intra-night variations in the optical band. A direct correla- 
tion with the emission at higher energies, which is not observed, 
might be washed out by the summation of several components. 
It is also possible that another component, e.g. emission from the 
accretion disc, is responsible for the rapidly varying optical flux. 

The observed optical variability on the time scale of several 
days (in the Watcher and ROTSE data) and the long-term varia- 
tion seen in the ROTSE flux might come from a different emis- 
sion region. For this component, emission regions of the size of 
the inner jet or large-scale jet are not excluded. The same regions 
might also account for the radio emission, which has been shown 
to be correlated with the optical flux on long time scales. 

4.2.3. Model 3 : Stratified jet scenario 

The third model considered here assumes a time-dependent strat- 
ified jet. Contrary to the previous models, this model has only 
one radiative component, the jet itself, the global SED resulting 
from the integration, along the jet, of the particle flow emission. 
A short description can be found in Appendix IB. 3 1 This model 
has alre ady been applied to describe the H .E.S.S. light curve of 
Flare 1 dBoutelier, Henri & Petrucci1l2008h . 

Fig. [13] shows the application of Model 3 to the data from 
Flare 2, with the principal model parameters summarized in 
Table lB.3l in Appendix lB.3l The model reproduces very well the 
SED of PKS 2155-304, as well as the simultaneous light curves 
in the X-ray and VHE range. It reproduces more accurately the 
variability in X-ray and VHE light curves than the previous mod- 



els, due to the freedom in rapidly varying the flux of particles 
from the injection and acceleration process, and it provides a 
full MWL prediction of the emission of the blazar. 

In the application presented here, the flux of injected parti- 
cles is kept roughly constant at the base of the jet, whereas the 
acceleration term varies. The flux of particles in the region where 
most of the VHE emission comes from is one order of magni- 
tude higher than the initial particle flux. This reveals that intense 
pair production has occurred between the base of the jet and the 
emission region, producing a strong amplification of initial per- 
turbations in the injection parameters, and hence the strong and 
fast variability observed at VHE. In the framework of Model 3, 
pair production appears to be a likely mechanism to interpret 
fast variability of PKS 2155-304. In situ re-acceleration of the 
electron-positron plasma and the additional emission from a first 
generation of pair-produced leptons in Model 3 provide the ob- 
served level of luminosity during Flare 2 with a smaller Doppler 
factor than in Model 2. 

The model partially reproduces the optical light curve over 
the night, although it does not fully match the amplitude of the 
observed flux variation. A smooth increase in the optical band 
is obtained as observed, in response to the rapid variation of the 
electron-positron flux. However it lasts for a shorter time inter- 
val than the one observed. This increase is due to the large size 
of the optical emitting region that plays the role of a low pass 
filter. The optical luminosity integrates the recent past history of 
the jet. No direct correlation between the optical and the high- 
energy flux is expected from the model, instead the optical flux 
responds to high-energy events in the source with a certain de- 
lay. Realistically, the optical component should contain contri- 
butions from earlier events that occurred before the onset of the 
observations during MJD 53946. Given a complete knowledge 
of the past variations of the high-energy flux, the inclusion of a 
longer history of high-energy events in this model might, at least 
in theory, lead to an exact match with the observed optical light 
curve. 
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Fig. 13. Time dependent SSC model (Model 3) for the night of MJD 53946. Left panel: SED of PKS 2155-304, including the band 
of the varying CHANDRA spectrum and the band of the varying H.E.S.S. spectrum. Nightly averaged spectra for CHANDRA and 
H.E.S.S. are indicated by red points. The optical flux detected by the Bronberg Observatory for this night is added as a red point; 
data from ATCA are presented as blue points. The thick line shows a snapshot of the modelled SED at a given state of the source, 
while the thin line shows the modelled spectra averaged over the H.E.S.S. observation period. Black data points indicate archival 
(non-simultaneous) data. Right panel: The points represent the H.E.S.S. (upper panel), CHANDRA (middle panel), and Bronberg 
(lower panel) light curves in units of 10" 10 erg cirT 2 s _1 . The solid lines are the light curves computed by the model in these three 
energy ranges. 



5. Discussion 

5.1. Modelling VHE flares and high-flux states 

Given the current understanding of the source PKS 2155-304, 
all three multi-zone SSC models presented in Section [4721 arrive 
at a good description of the VHE and X-ray emission during the 
night of Flare 2. This illustrates the present success of SSC mod- 
els for blazars, but also their inherent degeneracy. The model 
parameters are not unique solutions, but some general conclu- 
sions can still be reached. In all three scenarios, the dominant 
X-ray emitting zone appears larger than the VHE emitting zone. 
The X-ray emission can be attributed mostly to the jet, while the 
VHE emission comes from a denser emission region, identified 
as a blob in the jet or its inner core. The flares visible in the X- 
ray and VHE bands are ascribed to the dense emission region, 
with a size in the range of 10 15 to 10 16 cm and a magnetic field 
of a few 0.01 G. High Doppler factors (>30) are always needed 
and the high-energy tail of the particle distribution has to reach 
individual Lorentz factors of the order of 10 6 . 

A multi-zone scenario, similar to the ones described in 
Section|4] but on a longer time scale, might also account for the 
different characteristics observed during high and low states of 
the source. The synchrotron emission would be generally dom- 
inated by a not very dense emission region that does not con- 
tribute significantly to the SSC component. An additional emis- 
sion region with a denser and more energetic electron popula- 
tion would dominate the VHE emission while contributing lit- 
tle to the X-ray emission. This second zone would account for 
the steep correlation and spectral hardening during high states 
and flares. In a low state of the source, the small contribution 
from the second zone to the X-ray emission might become unde- 
tectable against the dominant emission from the first zone. This 



would naturally explain the uncorrelated behaviour or less steep 
correlations at lower flux states and the emergence of correla- 
tions during higher flux states. 

The three models presented in this work account for a large 
Compton dominance in the emission, yielding a VHE flux varia- 
tion in Flare 2 that by far supersedes the observed variation in the 
correlated X-ray flare. Clearly, an extreme situation could arise 
where the X-ray flux variation becomes too small to be detected 
against the emission from the larger zone and the VHE flare 
seems to be without any X-ray counterpart. Time-dependent 
multi-zone SSC scenarios can t hus account for the "o rphan 
flares" de tected in 1ES 1959+650 (|Krawczvnski et al. 120041) and 
Mrk 421 dBlazeiowski et al. 1120051). as previously s uggested by 
iKrawczvnski et al71 d2004l) and lLenain et al. I d2009l) . 

Complementary particle injection events can be added to 
Models 1 and 2 to reproduce more closely the structure of the 
light curves, which include at least one secondary peak during 
Flare 2. However, a description of the fastest observed variabil- 
ity during Flares 1 and 2 has not been attempted here. Variability 
of a few minutes puts strong constraints on the size of the emis- 
sion region and requires a very large bulk Lorentz factor of the 
emitting electrons, which is also neede d to allow the escape of 
VHE photons from the emission region dAharonian et al. | [2007: 
Begelman et al. 2008). Ascribing some of the VHE emission 
to external inverse Compto n processes might hel p to loosen 
some of these constraints (Begelman et al. 1 120081). Follo wing 
the "needle-in-jet" scenario (Ghisellini & Tavecchio 2008; i or 
the " minijet" scenarios (Giannio s et al.l 120091: iNalewaik o et al.l 
|201 it) , very small emission regions inside the jet, much smaller 
than the plasma blobs discussed in this work, could be respon- 
sible for very rapid variability. Such emission zones might be 
caused by magneto-centrifugal acceleration of electron beams 
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or magnetic reconnection in a Poynting flux dominated jet and 
would be moving rapidly to wards the observer. An alterna- 
tive e xplanation is given by iNeronov. Semikoz & Sibirvakov I 
(120081) . who assume a radiatively inefficient accretion flow and 
a relatively small mass of the central black hole to account for 
rapid variability. 

5.2. Low frequency emission and long-term evolution 

Additional constraints from MWL data should help to resolve 
the degeneracy of the models at high energies. However, despite 
their success at high energies, none of the above scenarios is able 
to provide a detailed description of the light curves at low fre- 
quencies. Qualitatively, Models 2 and 3 propose to interpret the 
optical emission as a low energy counterpart of the activity at 
high energies, with a dominant contribution from the extended 
synchrotron jet in Model 2, and integration of the recent past 
VHE events in Model 3. However, they do not reproduce quan- 
titatively the short-term optical light curve during Flare 2 or the 
long-term optical and radio evolution. This emphasizes the cur- 
rent difficulty to relate the very fast and sporadic events seen at 
VHE in blazars to long-term properties traditionally explored in 
the radio and optical ranges. Various clues show that they are 
partially linked, but how and to which extent remains unclear. 

In this regard, the radio and optical light curves gathered over 
years during and after the 2006 MWL campaign are quite inter- 
esting. The radio data might indicate a si milar situation to th e 
one seen in 2008 in the radio-galaxy M87 dAcciari et al. 1 120091) . 
In both cases, a rapid variation at the highest energies occurred 
during or shortly before an increase in the radio flux over several 
months. From the available data it seems that the radio flux from 
PKS 2155-304 was at a relatively high, but stable level when the 
VHE flares occurred, and started to increase about one month 
afterwards. 

In the two-zone scenarios of Models 1 and 2, part of the ra- 
dio and optical emission could come from the same lepton pop- 
ulation that causes the VHE y-ray and X-ray flares. The blob 
could be initially opaque to radio waves due t o synchrotron 
self-a bsorption, as has been observed in Mrk421 dCharlot et al. I 
2006). The expansion of the blob permits the emission of radio 
waves with a certain delay. The longer cooling times at lower 
energies lead to a longer time scale for radio and optical emis- 
sion compared to the very rapid high-energy flares. Moreover, 
the variability is diluted in the radio light curve by emission from 
the inner and large-scale jets, as well as from the extended radio 
lobes. Diffusion of energetic particles from the blob into an ex- 
ternal component could lead to a delayed response to the VHE 
flares in the low-energy bands. 

In the stratified jet of Model 3, the emission observed 
at lower energies represents the sum of different high-energy 
events, which are delayed and spread out in time due to the cool- 
ing, propagation and re-acceleration of leptons in the jet. As in 
the previous case, the plasma would be optically thick in the be- 
ginning and become transparent to radio synchrotron emission 
with adiabatic cooling. In this way, high activity at high ener- 
gies and on short time scales directly contributes to an averaged, 
long-term signature at low energies, while in the X-ray and VHE 
y-ray band, variations in the injected flux close to the base of the 
jet become visible as individual flares. 

For the observations of the VHE flare from M87 in 2008 
no simultaneous optical data have been published, thus it is not 
known if the optical flux showed a behaviour similar to the ra- 
dio flux. However, in the case of PKS 2155-304, the Watcher 
and ROTSE data show the average optical flux to be in a high 



state and on a long-term increase when the VHE flares occurred. 
Variations on the scale of a few nights and on intra-night scales 
are superimposed on this increase. Similar to the radio flux, the 
optical flux could represent a smoothed response to previous 
high-energy activity. 

The simulation of a much longer history of particle injection 
in the jet of Model 3, which is beyond the scope of the work 
presented here, would be needed to modify the optical and radio 
light curves and might account for the overall evolution in these 
wavebands. The sparse sample of VHE data does not permit to 
determine if the long-term increase in the optical range, already 
seen before the 2006 VHE flares, was triggered by an earlier 
state of high activity at high energies. Differences between the 
long-term radio and optical light curves might be explained by 
longer delays in the radio band from synchrotron self absorption 
and faster cooling times in the optical band. One additional dif- 
ficulty is that the short-term variations in the optical flux do not 
seem to be directly connected with variations at high energies. 
They could arise from the summation of a direct synchrotron 
component and a smoothed response to the higher energy bands. 
Optical polarization monitoring would be very useful in probin g 
these different components (e.g. lBarres de Almeida et a"fl l2010). 

5.3. Hadronic models 

An interesting alternative to leptonic scenarios as considered 
in this paper is given by hadronic models, where a population 
of relativistic hadrons, rather than electrons and positrons, is 
responsible for the dominant contribution to the high-energy 
emission. However, detailed time-dependent modelling with 
hadrons is difficult to achieve due to the higher complexity of 
the hadronic interactions and the large number of free param- 
eters in those models. Moreover, due to the low efficiency of 
the hadronic emission processes, such scenarios seem generally 
less adapted t o describe the high-energy emission from blazars 
(ISikoralboiOl) . 

An application of a hadronic scenario to the SED observed 
during Flare 2 would need to account for the close correla- 
tion between the X-ray and VHE flux, the dominance of the 
VHE emission, the (very) rapid variability time scale, and the 
low-frequency counterparts. Especially the reproduction of rapid 
variability presents a problem for hadronic models, since accel- 
eration and cooling time scales for hadrons are much longer than 
for electrons and positrons. A possible solution to this problem 
might be given by models where rapid variability is achieved 
through the inter action of baryonic m atter, e.g. in the form of 
stellar envelopes dBarkov et al. 11201 Oh . with the blazar jet. 

6. Conclusions 

A clear difference in the temporal evolution of the SED has 
been put in evidence between high- and low-activity states of 
PKS 2155-304. The analysis of correlation on various time 
scales of the flux evolution at different frequencies shows that 
the steep correlation between X-ray and VHE fluxes already seen 
during Flare 2 (MJD 53946) is also found in the evolution of the 
average fluxes over several nights with high VHE flux states. 
Spectral hardening with an increased flux, already found in the 
H.E.S.S. data, is also seen in the SWIFT XRT data, but not in 
the RXTE data, owing most likely to the harder energy band ob- 
served with RXTE. 

When the entire available MWL data of the 2006 campaign 
are considered, there is no compelling evidence of a correla- 
tion of the optical flux with the high state in the VHE range on 
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time scales of hours or days. However, the dramatic VHE flaring 
of PKS 2155-304 in 2006 occurred at a time when the source 
was already seen in a high-activity state at optical and radio fre- 
quencies. Moreover, the optical and radio luminosity further in- 
creased after the peak in VHE activity, although on a longer time 
scale. To the best of our knowledge, this is the first evidence of a 
relation between the long-term flux evolution at low frequencies 
(optical and radio) and VHE events in a blazar. This emphasizes 
the present need of global modelling of such sources to further 
explore the still missing link between their low and high-energy 
behaviour. 

The SEDs of the non-flaring nights of the MWL campaign 
in July 2006 are well reproduced by a stationary one-zone SSC 
model with only relatively small variations of the parameters that 
describe the 2003 low state. Time dependent multi-zone SSC 
scenarios arrive at good descriptions of the spectral and flux evo- 
lution in the X-ray and VHE bands during Flare 2. They are able 
to explain the observed dominance of the IC component in the 
high-energy bands and the steep correlation between VHE and 
X-ray fluxes with parameters in the generally accepted range 
for blazars and identify different emission zones. In general, it 
can be concluded that the SSC framework seems well adapted 
to describe high-energy data from blazars, even in the case of 
this very rich data set. However, these models still have diffi- 
culties to quantitatively reproduce the long-term low-frequency 
light curves. The more refined models allow in principle for a 
delayed reaction to the VHE flares in the optical and radio band 
as observed, but additional emission regions or mechanisms and 
further modelling of the long-term development of the source 
would be required to complete the picture and include the whole 
optical and radio emission in the simulations. Nevertheless, it 
remains to be seen if yet more refined multi-zone models will 
ultimately lead to a realistic description of the nature of the emis- 
sion, or whether the increasing number of free parameters pre- 
vent such conclusions, given the finiteness of observable data. 

Independent of the model details, however, a distinction 
of two physically different active VHE states of the blazar 
PKS 2155-304 is emerging, when considering the present data 
set together with previously published data on PKS 2155-304. 
In the low state, generally no correlation is seen between the 
X-ray and VHE flux, whereas the fluxes observed during high 
states show a steep correlation. Spectral variability is observed 
in the VHE band for low states, but it is different from the spec- 
tral hardening seen for high states. The different behaviour dur- 
ing low and high states might be explained within the two-zone 
scenarios considered in this paper. In this case, the synchrotron 
emission during low states would be dominated by a not very 
dense emission region, while a denser and more energetic elec- 
tron population would account for the emerging correlation and 
spectral hardening during high states. 
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Appendix A: Tables from the H.E.S.S., RXTE and SWIFT XRT spectral analyses 
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Table A.l. Average H.E.S.S. spectral parameters during the nights where X-ray data are available. The start time of observations To 
(in MJD), the flux normalization <po at 1 TeV (in 1CT 11 cirT 2 s -1 TeV" 1 ), the photon index F for a power law, the^ 2 over numbers 
of degree of freedom (ndof) and the live time of the selected runs (in ks) are liste d. * The entry with t he s tart time 53951.03 is an 
average over three nights (MJD 53951 to 53953). The quiescent flux is taken from Abramowski et al. I d2010l) . 
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Table A.2. RXTE spectral parameters from the different pointings of the nights during the H.E.S.S. campaign. The start time of 
observations To (in MJD), the extrapolated flux normalization <f>a at 1 keV (in crrT 2 s _1 keV~'), the photon index T for a power law, 
the^- 2 and ndof of the fit, and the live time of each pointing (in ks) are listed. 
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Table A.3. SWIFT-XKT spectra of the nights during and after the H.E.S.S. campaign. The start time of observations To (in MJD), 
the flux normalization <pa at 1 keV (in cirT 2 s keV -1 ), the photon indices below (Tsoft) an d above (Th ar j) a spectral break (Eb rea £ in 
keV), the x 2 over numbers of degree of freedom (ndof) and the live time of each pointing (in ks) are listed. Where no break energy 
is indicated, the third column provides the photon index over the whole spectrum. For the pointings on MJD 53945 and MJD 53952, 
the results of fits with a power law and with a broken power law with fixed break energy are shown. For the pointing on MJD 53953, 
a fit with a broken power law did not yield a satisfactory result, even when fixing the break energy. 
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Appendix B: A short description of the 
time-dependent models 

B.I. Model 1 

In Model 1 , a power-law distribution of relativistic electrons with 
index 2, a natural outcome of first order Fermi acceleration, is in- 
jected into the source volume, divided into several homogeneous 
cubic cells. The time development of these cells is simulated, 
including electron injection, radiative cooling, synchrotron and 
SSC emission, as well as absorption of light traversing the cells 
due to electron-positron pair production. Adiabatic expansion of 
the emission region is not taken into account. 

The model does account for the light crossing time effect be- 
tween the cells and a comoving observer, which becomes signif- 
icant in the case of rapid variability time scales. The comoving 
observer is assumed to view the jet under an angle of 90°, corre- 
sponding to a maximum delay of the light travel time for cells at 
different locations in the jet. The observed synchrotron spectrum 
is a superposition of spectra from different parts of the emission 
region and thus from electron distributions that have had more 
or less time to cool. This superposition leads to a break in the 
observed spectrum in a natural way, contrary to the stationary 
model discussed in Section 14.11 where a break is introduced a 
priori in the electron energy distribution. 

A more detailed description o f the model is g iven by 
Chiaberg e&Ghiseliinil d 1999ft and iKatarzvnski et al. I £2008). 
Table IB. II shows the set of parameters that were used for a de- 
scription of the varying SED during Flare 2. 
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Table B.l. Principal parameters used in Model 1 for a descrip- 
ti on of Flare 2. A definitio n of the parameters can be found 
in IKatarzvnski et al. I d2008l) . 



B.2. Model 2 

In this model, the high-energy emission comes from relativistic 
electrons and positrons in an extended, inhomogeneous jet and 
in a dense plasma blob travelling along the jet axis. 

The jet is at the origin of the long-term component, mod- 
elled as a steady radiation, that dominates the X-ray synchrotron 
emission. It has a global paraboloidal shape and is approximated 
by a continuous outflow of homogeneous cylindrical slices from 
the central engine, which travel with constant velocity and ex- 
pand adiabatically. The magnetic field intensity decreases along 
the jet. Adiabatic and synchrotron cooling and particle escape 



are accounted for, whereas the weak radiation field inside the 
extended jet makes IC cooling negligible. 

The inner plasma blob is a denser and more energetic zone, 
which is at the origin of the VHE flux and of the rapid variabil- 
ity through its SSC emission. The blob is built of homogeneous 
slices, in the same way as the jet. In the blob, the higher radiation 
field energy density implies a significant IC cooling, in addition 
to the cooling processes already considered for the particles in 
the jet. 

Radiative transfer along the jet is treated under the assump- 
tion of a small viewing angle in the comoving frame, thus the 
light crossing time effect described in Model 1 is negligible here. 
In addition to the radiation transfer inside each emission region, 
Model 2 also takes into account the absorption of synchrotron 
radiation from the blob by jet particles and vice versa, as well 
as the absorption of IC photons emitted from the blob through 
electron-positron pair production in the jet. A particle distribu- 
tion is injected at the base of the jet following a power-law with 
initial spectral index n ]eX in the jet and «hiob in the blob. 

Synchrotron emission by an additional component, identified 
as a large kpc-scale jet, dominates the radio and optical emission 
and completes the picture of the observed MWL emission from 
the blazar. Emission from this large-scale jet is treated here as 
independent from the other components, i.e. the inner jet and 
plasma blob. 

A more detailed descripti on o f this model i s given by 
IKatarzvnski. Sol&KusI (120031) and iLenainl d2009ft . Table lB~2l 
shows the set of parameters that were used for a description of 
the varying SED during Flare 2. 
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Table B.2. Principal parameters used in Model 2 for a descrip- 
tion of the high-energy emission of Fl are 2. A detailed descrip- 
tion o f the parameters can be found in Kat arzvnski. Sol & Kus I 
d2003l) . The parameters used for the large-scale jet that domi- 
nates the radio and optical emission are the same as given in 
Tablefflin Sectiongl] 



B.3. Model 3 

The third model fits well into th e t wo-flow fr amework origi- 
nally proposed bvlPelletierl (1 19851) and lSol et all d!989ft (see also 
iTsinganos & Bogovalov 2002 j or t he "spine-in-jet" model devel- 
opped by IChiaberge et al. 1 12000) where a non relativistic but 
powerful MHD jet launched by the accretion disk surrounds a 
highly relativistic plasma of electron-positron pairs propagating 
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along its axis. The MHD jet plays the role of a collimator and 
energy reservoir for the pair plasma, which is responsible for the 
observed broad-band emission. 

The jet has a paraboloidal shape and its magnetic field de- 
creases along the axis following a power law. A continuous ac- 
celeration of the plasma is assumed, leading to an increase of the 
bulk Doppler factor up to an asymptotic value dbcc, over an ac- 
celeration region starting at the base of the jet. The pair plasma 
is continuously being injected at the base of the jet, and its evo- 
lution is computed as it propagates inside the MHD structure. 
The plasma is re-accelerated along the jet to compensate for syn- 
chrotron and IC cooling, which are responsible for the observed 
high-energy emission, following an SSC scenario. Particle ac- 
celeration is supposed to be provided by the interaction of the 
plasma with the magnetic turbulence carried by the surround- 
ing MHD structure and is described as a power law of index f : 
Qdcc(z) °c 6o(0 {z/Zq)^. Hence, the particle energy distribution 
is a relativistic Maxwellian (or "pile-up") distribution, which is 
the natural outcome of second order Fermi acceleration. This is 
expected for e xample in the case of particle-wave interaction or 
reconnection dSchlickeiserfll984ll985HHenri & Pelletid[l99Tl 
ISauge & Henrill2004tlBoutelier. Henri & Petruccill2008l) . 

The evolution of the observed flux is reproduced as a conse- 
quence of time dependent injection parameters. Hence, by vary- 
ing the acceleration term at the base of the jet and the flux of 
injected particles, a time dependent SED and light curves in the 
different energy bands can be reproduced. Given the narrow en- 
ergy distribution and the effects of cooling of the plasma along 
the jet, the resulting SED is a superposition of contributions from 
different regions in the jet, where emissivity peaks at lower en- 
ergies for particles that have travelled farther. The high-energy 
peaks are dominated by emission close to the base of the jet. 

Contrary to the previously described two models, pair pro- 
duction is not only calculated to account for losses of VHE y- 
rays, but the first generation of electron-positron pairs is added 
to the particle flux and re-accelerated with the injected plasma. 
This leads to an enhancement of initial variations of the particle 
density, in particular during a flare. 

A more detailed description of this model is give n 
by iBoutelier, Henri & Petruccil d2008l) and iBoutelierl (120091) . 
Table IB. 31 shows the set of parameters that were used for a de- 
scription of the varying SED during Flare 2. 
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Table B.3. Principal parameters used in Model 3 for the de- 
scription of Flare 2. A definition of these parameters and 
details on constraints from obs erva tion can be fou nd in 
IBoutelier. Henri & Petruccil (T2008) and iBoutelier I (120091) . 
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